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Abstract 

Impairment of the tightly regulated ossification process leads to a wide range of skeletal dysplasias and deciphering their 
molecular bases has contributed to the understanding of this complex process. Here, we report a homozygous mutation In the 
mitochondria-associated granulocyte macrophage colony stimulating factor-signaling gene (MAGMAS) in a novel and severe 
spondylodysplastic dysplasia. MAGIVIAS, also referred to as PAM16 (presequence translocase-associated motor 16), Is a 
mitochondria-associated protein involved In preprotein translocation Into the matrix. We show that MAGMAS Is specifically 
expressed In trabecular bone and cartilage at early developmental stages and that the mutation leads to an instability of the 
protein. We further demonstrate that the mutation described here confers to yeast strains a temperature-sensitive phenotype, 
Impairs the import of mitochondrial matrix pre-proteins and induces cell death. The finding of deleterious MAGMAS mutations 
In an early lethal skeletal dysplasia supports a key role for this mitochondrial protein In the ossification process. 
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Introduction 

During embryonic development and postnatal growth period, 
most longitudinal bone growth occurs via endochondral ossification, 
a process in which chondrocytes produce a cartilage template that 
will ultimately be replaced by bone [1]. Disturbance of this tightly 
regulated process results in skeletal dysplasias (SD), which are an 
extremely diverse and complex group of rare genetic diseases [2] . 
Causative mutations have now been identified in over 230 diflFereiit 
genes in more than 400 unique skeletal pheiiotypes [3,4] . However, 
the genetic basis of over 100 different entities remains to be 
determined. Among SD, spondylodysplastic dysplasia is a generic 
descriptive term for a heterogeneous group of SD characterized by 
severe vertebral abnormalities and distinctive skeletal or extra- 
skeletal features [3,4] . The outcome is variable from early lethality 
to long-term survival. The last international classification recognizes 
at least five distinct entities within this group of rare disorders, i) 
achondrogenesis type lA [MIM 200600] due to mutations in 
TRIPll [MIM 604505], which encodes the thyroid hormone 
receptor interactor 11, ii) Schneckenbecken dysplasia [MIM 
296250] due to mutations in SLC35D1 [MIM 610804], encoding 
the solute carrier family 35 member Dl, also known as UDP- 
Galactose transporter-related 7 iii) spondylometaphyseal dysplasia 
(SMD) Sedaghatian type [MIM 250220], a heterogeneous condi- 



tion that is due in some cases, to mutations in SBDS [MIM 607444] 
involved in the processing of ribosomal RNA, iv) fibrochondrogen- 
esis [MIM 228520] , due in some cases to COLl lAl mutations and v) 
opsismodysplasia [MIM 258480], due to mutations in IKPPLl 
[MIM 600829] that encodes the inositol polyphosphate phospha- 
tase like- 1 [5] . The wide range of proteins involved in this group of 
SD points out the complexity of the ossification process. 

Although different entities within the group of spondylodysplas- 
tic dysplasia are well defined, a few cases remain unclassified. 
Recentiy, Megarbane et al. have reported on two unrelated 
consanguineous Lebanese families with four affected cases 
presenting a novel type of early lethal spondylometaphyseal 
dysplasia [6,7]. To elucidate the molecular basis of this entity, we 
performed exome sequencing in these two distinct families and 
identified a homozygous mutation in MAGMAS (Mitochondria- 
associated granulocyte macrophage colony stimulating factor- 
signaling molecule) in patients from both families. 

Results 

Patients 

Two unrelated consanguineous Lebanese families with four 
affected cases presenting a novel type of early lethal spondylodys- 
plastic dysplasia, recently reported by Megarbane et al. [6,7], were 
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Author Summary 

Skeletal dysplasias (SD) refer to a complex group of rare 
genetic disorders affecting the growth and development 
of the skeleton. The identification of the molecular basis of 
a considerable number of SD has greatly expanded our 
knowledge of the ossification process. Among SD, 
spondylodysplastic dysplasia is a generic term describing 
different conditions characterized by severe vertebral 
abnormalities and distinct by additional specific features. 
Several entities within this group are well defined. 
However, a few cases remain unclassified, of which a 
novel autosomal recessive spondylometaphyseal dysplasia 
recently reported by Megarbane et al. in two Lebanese 
families. Here, we identified MAGMAS as a candidate gene 
responsible for this severe SD. MAGMAS, also referred to as 
PAM16, is a mitochondria-associated protein, involved in 
pre-proteins import into mitochondria and essential for 
cell growth and development. We demonstrated that 
MAGMAS is expressed in bone and cartilage in early 
developmental stages underlining its specific role in 
skeletogenesis. We also give strong evidence of the 
deleterious effect of the identified mutation on the in- 
vivo activity of Magmas and the viability of yeast strains. 
Reporting deleterious MAGMAS mutation in a SD supports 
a key and specific role for this mitochondrial protein in 
ossification. 



included in tliis study. The main clinical features were pre- and 
postnatal growth retardation, developmental delay, dysmorphic 
features, prominent abdomen, narrow thorax with short ribs and 
short limbs. Skeletal features included severe platyspondyly at 
birth that improved somewhat over time, square iliac bones, 
horizontal acetabulae, short long bones with abnormal modeling, 
widening of the distal femoral metaphyses and delayed epiphyseal 
ossification (Figure 1). In the first reported family (Family Fl, 
Figure 2) [6], both affected sibs died at an early age from 
respiratory insufficiency (Fl-IV. 1 at age 9 months and Fl-FV.S at 
age 2 years) while heart failure was the cause of the early death, at 
age 2 years, in patients (F2-IV.2 and F2-IV.3) of the second family 
(Family F2: individuals IV.2 and IV.3) [7] (Figure 2). 

Molecular findings 

Exome sequencing has been performed in Family 2, in three 
individuals: the proband (F2-IV.3) and her two unaffected parents 
(F2-in.l and F2-III.2) (Figure 2). Considering the recessive 
inheritance of the disease in the pedigree and consanguinity 
(F = 1 / 1 6), we assumed identity by Descent [8] . Using an in-house 
bioinformatics tool, we realized identity by descent filtering of 
exome data and identified 1037 homozygous variants in the 
proband, which were also heterozygous in both parents. Of these 
1037 variants, 916 were exonic and/or splice site mutations, of 
which 434 were either non-synonymous and nonsense single 
nucleotide variants, or small insertions/ deletions leading to a shift 
in the reading frame. 

In order to further refine the list of candidate nucleotide 
variations, additional filtering was performed on what we called 
the "dbSNP13700" database (i.e. dbSNP137, from which we have 
removed mutations that are reported as pathogenic alleles, as well 
as all referenced variants with frequencies lower than 1%, or 
without any frequency information). Additional filtering was also 
performed using our in-house exome database that includes all 
homozygous variants detected by exome sequencing of 14 
Lebanese unaffected individuals. We ended up with a fist of 16 



candidate homozygous by descent variants in the proband F2- 
IV. 3. 

In family 1 , exome sequencing has been performed only in the 
proband (F1-IV.3, Figure 2). Using the same analysis pipeline and 
after all filtering steps, we were able to identify 59 homozygous 
variants in the genes coding regions of this patient. 

In total, patients F1-IV.3 and F2-IV.3 shared 13872 homozy- 
gous variants, of which only 3 exonic variations were both 
homozygous by descent in the proband from family 2 (F2-IV.3), 
and homozygous in the proband from family 1 (Fl-IV. 3), after 
applying the different filters (Table SI). Of those three variants, 
only one missense mutation (NM_016069: c.226A>G; p.As- 
n76Asp) in the PAM16 gene (presequence translocase-associated 
motor 16), that is more commonly referred to as MAGMAS 
(Mitochondria-associated granulocyte macrophage CSF-signaling 
molecule), was predicted to have pathogenic effects using different 
bioinformatics prediction softwares, such as: Mutation Taster 
(disease causing variation with a probability value of 0.99999957), 
SIFT (damaging effect, score 0.006) [9], PROVEAN (deleterious 
effect, score —4.021) [10] and PolyPhen-2 (possibly damaging 
effect, score 0.728 HumDiv). By Sanger sequencing, we confirmed 
the segregation of the c.226A>G transition in MAGMAS vnlh. the 
disease in both families (Figure 3A): the mutation was homozygous 
in the patients, heterozygous in the parents and in the unaffected 
sibs in both families. We also demonstrated that the variation was 
absent in a set of 550 Lebanese control chromosomes. Protein 
sequence alignments showed that the Asparagine at codon 76 in 
MAGMAS is well conserved among mammals and vertebrates 
(Figure 3B). 

The two patients described here are not related. However, they 
display the same homozygous missense mutation, thus suggesting 
the existence of a founder mutation. In order to verify this 
hypothesis, we used both exome data and genotyping of STR 
markers. By analyzing exome data, we were able to identify 1 1 
regions of shared homozygosity on chromosome 16. However, 
only three regions displayed a size over 100 kb, of whom two were 
containing or contiguous to the MAGMAS\oa\i (Chrl6: 4015729- 
7102036 (~3Mb) and Chrl6: 8728764-8840025 (-1.1 Mb)). 
Also, by using PLINK [11], we were able to reconstruct 
haplotypes from exome data of both patients, and here again, 
we identified a 3.6 Mb shared homozygous region located 
between SNP markers rs37301 19 and rs4616299 on chromosome 
16pl3.3 (Chrl6:4,057,603-7,657,432). Genotyping of microsatel- 
lite markers allowed to refine this region and define a minimal 
ancestral common homozygous haplotype spanning 1.9 Mb 
between markers D16S758 and D16S243 containing MAGMAS 
(Figure SI). Altogether, our results demonstrate that the c.226A> 
G transition in iM/lGA£45(p.Asn76Asp) is likely to be pathogenic in 
both families. MAGMAS \i located on chromosome 16pl3.3 and is 
composed of five exons. It encodes MAGMAS, an essential 
component of the pre-sequence translocase-associated protein 
import motor (PAM) that regulates pre-protein translocation into 
the mitochondrial matrix [12]. 

Gene expression analysis 

A naturally occurring read-through transcription between the 
neighboring C0R07 (coronin 7) and AIAGAIAS genes on 
chromosome 16 is stated in the Genome UCSC browser 
(NM_00 120 1479). To determine if the fusion protein is involved 
in the pathogenesis of the severe skeletal manifestations observed 
in the affected cases, we checked the expression of C0R07, 
MAGMAS and C0R07-MAGMAS fusion transcript in human 
cultured control cells (fetal chondrocytes, osteoblasts and skin 
fibroblasts) by Reverse Transcriptase PGR. C0R07, MAGMAS 
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Figure 1 . Radiological features of the patients F2-IV.3 and F1 -IV.3. Radiographs of patient F2-IV.3 at 9 months (A, C, E) and F1 -IV.3 at birth 
(B, F) and at 3 months (D) show platyspondyly, square iliac bones, and delayed epiphyseal ossification. 
doi:1 0.1 371 /journal.pgen.1 00431 l.gOOl 



PLOS Genetics | www.plosgenetics.org 



3 



May 2014 | Volume 10 | Issue 5 | e1004311 



MAGMAS and Endochondral Ossification 



F1 



F2 



a 



IV 



a 



o 



IV 



O 




6, M 0 



Figure 2. Pedigrees of the families F1 and F2 included in this study. Blackened symbols indicate affected individuals. 
doi:1 0.1 371/journal.pgen.1 00431 1 .g002 



and C0R07-MAGMAS were all expressed in skin fibroblasts. 
MAGMAS expression in bone and cartilage was confirmed. On die 
other hand, neither C0R07 nor C0R07-MAGMAS was expressed 
in significant amounts in chondrocytes and osteoblasts (Figure 4A). 
This result was confirmed by quantitative real-time RT-PCR 
(Figure 4B): The ratio of MAGMAS in chondrocytes versus in 
fibroblasts is 0.82 whereas the ratios of C0R07 and C0R07- 
MAGMAS in chondrocytes versus in fibroblasts are 0.24 and 0.12, 
respectively. 

Immunolocalization of MAGMAS protein in wild type 
growth plates 

To further confirm the specific involvement of MAGMAS in 
skeletogenesis, we performed immunohistochemical staining of 
fixed femurs isolated from wild type (WT) mice at different 
developmental stages from El 6.5 to 2 weeks of age. Type X 
collagen staining was also performed in order to mark the 
terminally differentiated hypertrophic chondrocytes. MAGMAS 
protein was expressed in bone and cartilage. Expression of the 
protein was mostiy observed in the hypertrophic zone of the 
cartilage mainly at postnatal ages and in the primary spongiosa of 
trabecular bone. Its expression in proliferative chondrocytes and in 
the resting zone of the cartilage was not detected (Figure 5). 



medium (Figure 6A, 6B). To figure out if the Asn76Asp mutant 
confers a deleterious growth phenotype at non-permissive 
temperatures, tenfold serial dilutions oi pamlGA spores expressing 
WT or mutant MAGMAS were subsequently spotted onto 
fermentable (YPD) and non-fermentable (YPG) growth media 
and incubated at different temperatures for 3 days. As shown in 
Figure 6C, the mutant strain ^am7 6/4 -MAG^sj^yg^^p that presented 
a normal growth phenotype on both YPD and YPG media at 
28°C (301.I5°K), grew more slowly than paml6A-MAGy„ on 
YPD at 34°C (307.15°K) and 36°C (309.15°K). A more severe 
growth defect of /)ami6'.4-MAGAs„76As|) on YPG was evidenced by 
a very slow-growth at 34°C (307.1,5°K), and inviability at 36°C 
(309.15°K). To confirm this temperature-sensitive growth pheno- 
type, we secondly tested the ability oi MAGMAS j^^ygj^p to rescue 
the inviability otthe paml6A strain by plasmid shuffling using the 
haploid YP}i499paml 6A strain transformed with a plasmid 
encoding the yeast WT RAMI 6 iyPAMlG). Unlike WT MAGMAS 
expressing strains, the mutant /)flm76'zl-MAGA.,ii76As|> exhibited an 
obvious slow-growth at 34°C and was unable to grow at 36°G 
(309.15°K) on YPD (Figure 6D). On die otiier hand, paml6A- 
PAMi ff-MAGAsnTSAsp c^Us Carrying a WT copy of yPAMlG in 
addition to the mutant MAGMAS expressing plasmid, showed a 
normal growth phenotype at aU temperatures. 



Functional complementation 

Human and murine forms of MAGMAS present a very high 
degree of homolog}'. Of the 125 amino acids comprising the 
mature protein, 120 are identical [13]. The evolutionary 
conservation oi MAGMAS from yeast to mammals [14] and the 
lethal phenotype of its deletion (paml6A) in Saccharomyces cerevisae 
[15] and Drosophih [16] imply a fundamental role of this protein in 
cell growth and development. Previous studies have shown that 
human MAGMAS exhibits a complete growth complementation of 
deleted /)am7 6' yeast cells at all temperatures [17]. 

In order to obtain in vivo evidence of the deleterious nature of 
the identified mutation, we performed a yeast complementation 
assay using genetic approaches as described previously by Sinha et 
al. [17]. We first found tiiat at 28°C (301.15°K), bodi WT and 
mutant MAGMAS genes were able to complement the growth 
defect of the pamlSA cells since spores expressing iM/4GAfc46'fyr and 
MAGMASAm76Asp were viable when spotted onto G418 selective 



Protein import and MAGMAS expression in yeast extracts 

MAGMAS is an essential component of PAM, a peripheral 

mitochondrial inner membrane complex, that is crucial for 
preprotein translocation into the mitochondrial matrix [12,18]. 
Preproteins targeted to the matrix cross the outer mitochondrial 
membrane through the TOM complex then pass through the 
inner membrane via the TIM23^*^™''^ complex which is 
composed of the TIM23 translocase associated with the PAM 
complex [19] (Figure S2). To assess the transport function of the 
mutant MAGMASAsnveAsp, we evaluated the protein translocation 
activity by following the mitochondrial import of the HspGO 
protein. Non-processed mitochondrial HspGO was accumulated in 
significant amounts in the /)fflm76'/l-MAGAsn76Asp strain when the 
mutant phenotype was induced (37°C, 310.15°K), but was 
undetectable in WT!" paml 6A-PAA4} 6 and paml 6A -MAGMAS 
expressing cells (Figure 7A). We then evaluated the expression 
levels of MAGMAS protein in the difierent strains, by performing 
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Figure 3. The missense MAGMAS mutation c.226A>G (p.Asn76Asp). (A) Sequencing chromatograms showing the segregation of the 
c.226A>G transition in MAGMAS with the disease in both families F1 and F2. n: heterozygous peak. (B) Multiple alignments between human 
MAGMAS protein and several orthologs showing that the Asparagine at codon 76 in MAGMAS is well conserved among mammals and vertebrates. 
doi:1 0.1 371 /journal.pgen.1 00431 l.gOOS 



a Western blot using Human anti-MAGMAS antibodies. Immu- 
noblotting detected a lower level of mutant MAGMAS when the 
strain was shifted at the non permissive temperature, compared to 
the amount of MAGMAS^t ™ the same condition (Figure 7B). 

Mitochondrial morphology and peroxisomal biogenesis 
defects in the MAGMASAsn76Asp expressing strain 

Previous studies in mammalian cells, yeast and drosophila have 
shown the crucial role oi MAGMAS and its respective homologs in 
mitochondria biogenesis [12,17,20]. Altered mitochondria mor- 
phology observed in /(flm/ff yeast mutants and in 5//)-knowdown 
drosophila cells led us to study the aspect of these organelles in our 
mutant strains [20,21]. Microscopic examination of mitochondria 
revealed an unhealthy feature in the mutant strain consisting in 
punctuate and highly fragmented mitochondria that differed from 
the normal reticulated mitochondrial network observed in WT 



MAGMAS expressing cells. This unusual aspect was more 
obviously seen in the mutant strains shifted at the non permissive 
temperature (Figure 8A). 

Furthermore, the recent finding of impaired peroxisomes 
biogenesis in paml6 yeast mutants [21] prompted us to check 
the morphology of these organelles in our mutant strains. A 
remarkable aspect (vacuole labeling), indicative of pexophagic 
degradation, was only observed in pamWA-MAGp^^^-^-if^p^^ cells 
grown at the non permissive temperature compared to pamlSA- 
MAGwT cells (Figure 8B). 

Discussion 

Here, we report the identification of a homozygous missense 
AiiGAMiS mutation, (p.Asn76Asp), in patients from two unrelated 
Lebanese families, affected with a rare lethal spondylometaphyseal 
dysplasia recently described by Megarbane et al. [6,7]. By 
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Figure 4. Transcript expression analysis of MAGMAS, C0R07 and C0R07-MAGMAS. (A) Expression analysis oi MAGMAS, COfiOZand C0R07- 
MAGMAS by RT-PCR in control fibroblasts, chondrocytes and osteoblasts, showing that, contrary to MAGMAS, neither C0R07 nor C0R07-MAGMAS 
was significantly expressed in chondrocytes and osteoblasts. (B) Quantitative expression analysis of MAGMAS, C0R07 and C0R07-MAGMAS showing a 
very low level of expression of C0R07 and C0R07-MAGMAS in chondrocytes compared to that of MAGMAS in the same sample. 
doi:1 0.1 371 /journal.pgen.1 00431 1.g004 



reconstructing haplotypes from exome and STR genotyping data 
from patients F1-IV.3 and F2-IV.3, we were able to identify a 
minimal common ancestral homozygous haplotype around 
MAGMAS, spanning 1.9 Mb on chromosome 16pl3.3, suggesting 
the existence of a founder mutation in gene. MAGMAS, also 
referred to as PAA416, was first identified in 2001 byjubinsky et al. 
[13], as a gene encoding a novel mitochondria-associated protein 
involved in granulocyte macrophage colony stimulating factor 
(GM-CSF) signal transduction and expressed by nearly all 
mammalian cells. The observation of a AlAGAdAS mutation in a 
severe chondrodysplasia supports a specific role of MAGMAS 
protein in the endochondral ossification process. Its expression in 
the notochord and the somites has been reported at day 12.5 
embryo and in the cartilage of the developing ribs at day 14.5 
embryo [22]. However, its involvement in skeletogenesis has never 
been described to date. Our findings demonstrate that MAGMAS is 
expressed in trabecular bone and cartilage and more specifically 
by differentiated chondrocytes localized in the hypertophic zone 
and by osteoblasts at early developmental stages underlining its 
key role in skeletogenesis. On the other hand, the implication of 
the fusion protein C0R07-MAGMAS in the pathogenesis of the 
skeletal dysplasia has been ruled out, as the related transcript was 
not well expressed in human fetal chondrocytes and human 
osteoblasts. 

Besides the skeletal features, respiratory insufficiency (Fl-IV. 1 
and Fl-rV.3) or heart failure (F2-IV.2 and F2-IV.3) were 



observed, leading to early deaths of the patients [6,7]. Of note, 
an elevated MAGMAS level was detected in murine heart and lung 
in previous studies [22]. These features may also support a key role 
of MAGMAS in heart and lung. 

MAGMAS consists of two main functional domains: a N- 
terminal hydrophobic region including a membrane association 
domain (M) and a predicted mitochondria targeting region (T), 
and a C-terminal J-like domain (J) characterized by three ot-helical 
segments [12,23]. Importantly, the p.Asn76Asp mutation lies in 
the helix II of MAGMAS J-like domain. The latter domain forms, 
with the J-domain of DNAJC 1 9 (ortholog of yeast Pam 1 8 which is 
another essential component of PAM), a stable heterodimeric 
subcomplex [19,24,25] (Figure 82). Mutations in DJVAJC19 
preventing MAGMAS/DNAJC 19 dimerization have been shown 
to lead to dilated cardiomyopathy with ataxia syndrome (DCMA) 
[26,27]. Considering the localization of the p.Asn76Asp mutation 
in the J-like domain of MAGMAS, a similar phenotype was likely 
to be observed in our reported patients. Cardiomyopathy was 
observed exclusively in patients of the second family and ataxia 
was not present but all 4 cases died early before 2 years of age. On 
the other hand, the absence of DCMA in affected cases may also 
suggest a redundancy in the function of the subcomplex 
MAGMAS/DNAJC 19 in some tissues. 

We further demonstrated that the p.Asn76Asp mutation confers 
to yeast strains a temperature-sensitive phenotype characterized by 
a very slow-growth at 34°C (307.15°K) and inviability at 36°C 
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Figure 5. MAGMAS localization in control growth plates of mice. Distal femurs of wild-type mice at different developmental stages; (A) 16.5- 
day embryo; (B) Newborn and (C) 2 week-old; were fixed and stained with anti-MAGMAS antibodies. Sections were also stained with anti-CoX (marker 
of differentiated hypertrophic chondrocytes) and with HES. MAGMAS expression is reflected by the brown precipitate resulting from the peroxidase 
reaction. MAGMAS was detected mainly in hypertrophic chondrocytes (black arrows) and osteocytes of control growth plates at the different 
developmental stages. Its expression was also detected in some type X collagen -negative chondrocytes that are initiating the process of 
differentiation. 

doi:1 0.1 371 /journal.pgen.1 00431 l.gOOS 



(309.15°K), on both fermentable and non-fermentable growth 
media. Remarkably, the deleterious effect of the mutation was only 
observed at a homozygous state which is in correlation with the 
recessive mode of inheritance of the spondylody.splastic dysplasia 
reported in this study. 

Our findings of impaired import of mitochondrial-matrix 
proteins in the mutant strains shifted at non-permissive temper- 
ature confirmed the deleterious nature of the identified mutation. 
Likewise, it has been demonstrated that aa alterations in 
MAGMAS J-like domain resulted in temperature sensitivity and 
in vivo protein import defects in yeast cells [17,28]. The low 
expression level of MAGMAS detected in the mutant strain at the 
non permissive temperature indicates that the p.Asn76Asp 
mutation confers protein instability. Additionally, we observed 
unusual punctuated and highly fragmented mitochondria together 
with enhanced pexophagy in the mutant MAGMAS strain at the 
non permissive temperature, supporting a cell death induction in 
the yeast mutant strain. These data are consistent with the finding 
of Roy et al., in Drosophila [20] and confirm the crucial role of 
MAGMAS for cell integrity and viability. 

The association between a mitochondrial protein and a skeletal 
dysplasia is striking as skeletal manifestations are rarely seen in 
mitochondrial disorders. To date, the only severe skeletal dysplasia 
identified as having a possible mitochondrial basis is the anauxetic 
dysplasia [MIM 607095], which is allelic to two milder forms of 
SD, namely metaphyseal dysplasia Mc Kusick type [MIM 250250] 
and metaphyseal dysplasia without hypotrichosis [MIM 250460]). 
The gene responsible for these SD is RMRP [MIM 157660], which 
encodes the RNA component of the mitochondrial RNA 
processing endoribonuclease (MRP) complex [29-31]. Of note, 
RNase MRP is not exclusively found in the mitochondria but also 
in the nucleus. Previous studies have shown that mitochondrial 



and nucleolar RNase MRP have identical RNA components but 
distinct enzymatic activities and protein components [32]. The 
recent identification of mutations in POPl (Processing of precursor 
1) [MIM 602486], which encodes a protein component of the 
nuclear ribonucleoprotein complex in two siblings with anauxetic 
dysplasia [33] may support the involvement of the nuclear 
ribonucleoprotein complex rather than the mitochondrial complex 
in the etiology of these severe skeletal dysplasias. Similarly, besides 
its role in preprotein import into mitochondria. Magmas may also 
be involved in additional important mechanisms as suggested by 
Short et al. [21]. The detection of different human transcripts of 
MAGMAS, with some lacking the mitochondrial transit peptide 
(mTP) leader sequence [14] may support a novel extra mitochon- 
drial activity of this protein. 

In conclusion, this is the first report linking a mutation in 
MAGMAS v/ilh a human developmental disorder and supporting a 
key role for this mitochondrial protein in the ossification process. 

Materials and Methods 

Patients 

Informed consents for participation, sample collection and 
photographs publication were obtained from both families in 
compliance with the ethics guidelines of "ConseU de recherche de 
rUniversite Saint-Joseph" (number PTS5). Peripheral blood 
samples were collected from all available family members. 

DNA extraction and exome sequencing 

Genomic DNA was extracted from peripheral blood samples by 
standard salt-precipitation methods [34]. Exome sequencing was 
carried-out on 3 individuals from family 2 (F2-III.1, F2-III.2 and 
F2-IV.3) and the proband from family 1 (individual F1-IV.3) 
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Figure 6. In wVo growth analysis. (A-C) Diploid BY4743z(//p3/+ cells were transformed with (A) pRSAT 5-WT-MAGMAS and (B) pRS41 5-Asn76Asp- 
MAGMAS, and then subjected to sporulation followed by tetrad dissection analysis. Spores from a single ascus grown on YPD were spotted onto: 
G418 for the selection of spores carrying the deleted paml6 gene {pam 1 6A::kanMX6) and onto Leu drop-out media for the selection of MAGMAS 
positive spores. (C) Growth complementation analyses showing the temperature-sensitive growth phenotype of the mutant pam16A-MAG^s„76/isp in 
the BY4743 genetic background at 34 C (307.1 5"K) and 36"C (309.1 5"K) on both fermentative (YPD) and non-fermentative (YPG) media. (D) Hapioid 
YPH499 paml6A/pPAM16 cells transformed with pRS41 5-WT-MAGMAS and pRS41 5-Asn76Asp-MAGMAS were tested for their ability to rescue the 
inviability of the paml6A strain by plasmid shuffling at different temperatures. 
dol:l 0.1 371/journal.pgen.1 00431 1 .g006 



(Figure 2). Targeted exome sequencing, library preparation, 
capture and sequencing were performed by the French National 
Genotyping Center (CNG). Exomes were captured and enriched 
using the in solution Agilent SureSelect Human All Exon kit v3.() 
[35] and then sequenced on an lUumina HiSeq2000, using a 
paired-end 100-bp read sequencing protocol. Image analysis and 
base calling were performed using the Illumina Data Analysis 
Pipeline Software 1.5 with default parameters. Raw data were 
mapped to the current built of the human genome (hgl9) by using 
BWA 0.75 [36] . Variant calling was subsequently performed using 
GATK 2.5.2 [37] and annotation was done with ANNOVAR 
[38]. In order to predict the deleterious effect of the identified 
sequence variations, different bioinformatics tools were applied; 
such as MutationTaster (http://w\vw.mutationtaster.org/), SIFT 
(http://sift.bii.a-star.edu.sg/) [9], PROVEAN (http://provean. 
jcvi.org/index.php/) [10] and PolyPhen-2 (http://genetics.bwh. 
harvard.edu/ pph2/). 

Linkage analysis 

Searching for shared homozygous regions using exome data and 
reconstructing haplotypes were realized using the free software 



package PLINK (http://pngu.mgh.harvard.edu/purcell/pUnk/) 
[11]. 

Genotyping of SIR markers 

The following STR markers flanking the PAM16 gene and 
spanning a 8.7 Mb region on chromosome 16pl3.3, were amp- 
lified using fluorescently labelled primers: D16S521 (AFMAl 
39WG1), D16S3024 (AFMAl 34XC1), D16S3070 (AFMA33Y 
HI), D16S3027 (AFMAl 54WC9), D16S758 (CHLC.ATA24C09), 
D16S3084 (AFMB080YH5), D16S3072 (AFMB015WA9), D16S 
3134 (AFMA059WF9), D16S510 (AFM312VD5), D16S423 (AFM 
249YG5), D16S423 (AFM249YC5), AFM339XG1 (RH4035) and 
D16S3135 (AFMA059XC1). Marker heterozygosity rate ranged 
from 0.71 (D16S521 and DI6S51()) to 0.87 (D16S3082 and 
D16S3027). PGR amplifications were performed in a 30-|il 
reaction mix, containing 50 ng of genomic DNA, 0.5 units of 
Taq polymerase (Invitrogen), 1.5 mM MgCl2, 5 pmol of each 
primer, and 0.2 mM of each dNTP. After an initial dena- 
turation step of 5 min at 96°C, 35 cycles of amplification (94°C for 
30 sec; 55°C for 30 sec; 72°C for 30 sec) were performed, 
followed by a final elongation step of 5 min at 72°C. After 
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Figure 7. Preprotein translocation and iVIAGMAS expression analyses in yeast cells. paml6A-MAG^ and paml eA-MAG^snzsAsp strains in 
the YPH499 genetic context, expressing wild-type and mutant MAGMAS, accordingly were constructed by plasmid shuffling (see Materials and 
Methods). Strains were grown at 28 C (301.15'K) then shifted for 8 h (28800 s) at 37X (310.15 K) to induce the mutant phenotype. (A) Analysis of 
protein translocation in yeast cells: Non-processed precursor form and mature Hsp60 were detected by Western blotting using anti-Hsp60 antibodies. 
The non-processed Hsp60 precursor accumulated in the mutant strain after induction of the phenotype. Only the mature form of Hsp60 was 
detected in wild-type yP/\/M;6 and MAGMAS expressing yeast cells. PGK was used as a loading control and Porin as a mitochondrial loading control, p: 
precursor form of Hsp60, m: mature form of Hsp60. (B) Immunoblotting using Human anti-MAGMAS antibodies detected a lower level of MAGMAS 
protein in the mutant strain when shifted at non-permissive temperature, compared to the amount in MAGMASwt expressing strains. 
doi:1 0.1 371 /journal.pgen.1 00431 1.g007 



amplification, PGR products were subjected to denaturing 
electrophoresis on an ABI PRISM 3100 DNA Sequencer (Applied 
Biosystems). 

Capillary sequencing 

Genomic and cDNA sequences of AIAGAiAS (also known as 
PAM16) were obtained from UCSC Genomic Browser on Human 



February 2009 (NM_0 16069). Primers used for PGR amplification 
were designed using PrimerS software (http://frodo.wi.mit.edu) to 
amplify the region surrounding the mutation detected by exome 
sequencing in exon 4. PGR products were purified by exonuclease 
I/Shrimp Alakaline Phosphatase treatment (ExoSAP-IT; Fisher 
Scientific SAS, lUldrch, France) according to the manufacturer's 
instructions and both strands were sequenced using the Big Dye 
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Figure 8. Morphological defects of the yeast MAGMAS nwA^nX. strain. (A) For mitochondria visualization, YPH499-pam76z(-MAGv(n- and 
YPH499-pamy6zl-MAGAsn76Asp strains were transformed with pYX232-mtGFP plasmid; microscopy filters used were 450/490 nm for excitation and 
500/550 nm for emission. Unlike the normal reticulated mitochondrial network observed in wild-type MAGMAS expressing cells, highly fragmented 
mitochondria were detected in the mutant strains grown at both permissive and non-permissive temperatures. (B) For peroxisomes visualization, BY- 
pam16A-MAG<j^ and BY-pam;6z/-MAG^_s„76^5p were transformed with the pUG34DsRed.SKL reporter plasmid; filters were 542/582 nm for excitation 
and 604/644 nm for emission. A red halo (white arrow) indicative of pexophagic degradation was observed by mutant cells when grown at the non 
permissive temperature. 
doi:1 0.1 371/journal.pgen.1 00431 l.gOOS 



Terminator vl.l Cycle Sequencing Kit (Applied Biosystems). 
Sequence reactions were purified on Sephadex G50 (Amersham 
Pharmacia Biotech, Foster City, CA) and capillary electrophoresis 
was performed on Genetic Analyser 3100 (Applied Biosystems). 
Electropherograms were analyzed on the Sequence Analysis 
Software version 5.2 (Apphed Biosystems) and aligned with the 
WT MAGAIAS gene sequence using ChromasPro version 1.22 
(Technelysium, Queensland, Australia). 

Reverse Transcriptase PCR 

Total RNAs were extracted from human primary cultured 
control cells (fetal chondrocytes, osteoblasts and skin fibroblasts) 
via the RNeasy Mini Kit (Qjagen). Complementary DNA (cDNA) 
was synthesized using the Superscript II Reverse Transcriptase 
(Invitrogen Life Technologies, Carlsbad, CA, USA) according to 
the manufacturer's protocol. Sequences of MAGMAS 
(NM_016069), C0R07 (NM_001201472), C0R07-MAGMAS 
(NM_00 120 1479) and (NM_002046) were obtained from 

UCSC Genomic Browser on Human February 2009 and specific 
primers for cDNA amplification were designed through the 
PrimerS software as follows: 

MAGMAS-Forward: 5 '-TACCTGGCCCAGATCATTGT-3 ' , 
MAGMAS-Reverse: 5'-GCATCTGCCCTTTTTCTCTG-3'; C- 
OR07-Forward: 5'-GTCCTGGTGTACTGGGCATT-3', CO- 
R07-Reverse: 5'-AGAGTCAGGGTCCAGCAGAG-3'; COR- 
07-MAGMAS-Forward: 5'-CTCAGCGCAGTACCTGGAAG- 
3', C0R07-MAGMAS-Reverse: 5'-GCCACCCACGGATTT- 
ATCAT-3' and GAPDH-F: 5'-ATGTTCGTCATGGGTGTG- 
AA-3', GAPDH-R: 5'- TTCAGCTCAGGGATGACCTT-3'. 



Quantitative Real-Time RT-PCR 

Reactions were conducted in a 96-well plate with the ABI 7500 
Sequence Detection System (Applied Biosystems). Primers against 
MAGMAS, C0R07 and C0R07-AUGMAS were designed using 
Primer Express software (Applied Biosystems) (Figure S3). PCR 
was performed in a 20-|Xl reaction volume containing 1 0 |xl Power 
SYBR-GreenPCR Master Mix of SYBR Green PCR Master Mix 
buffer (Applied Biosystems, Foster City, CA, USA), 10 pmol 
forward and reverse primers and 5 ng of RNAse H-treated 
(Invitrogen Life Technologies) cDNAs obtained from reverse 
transcription of DNAse I-treated total RNA extracted from control 
fibroblasts and chondrocytes cDNA. The reaction cycling 
conditions were 95°C for 10 min, followed by 40 cycles of 95°C 
for 15 s and at 60°C for 1 min. Sequence Detection Software was 
used for exporting the threshold cycle data and further analyzing 
differences in threshold cycle values (ACt) between the test locus 
and the control locus. The expression values are normalized to 
that of GAPDH. For each gene, the expression ratio for a sample 
was calculated as the ratio between the average mRNA signal in 
the sample and the signal from fibroblasts 

Immunohistochemical staining 

Femurs isolated from WT mice (at diflFerent developmental stages 
including E 16.5-day, birth and 2 weeks of age) were frxed with 4% 
paraformaldehyde and embedded in paraffin. Femur sections were 
stained with Hematoxylin and eosin (H&E) using standard protocol 
for histological analysis or were subjected to immunohistochemical 
staining. For immunohistochemistry, sections were stained with Ab 
specific to MAGMAS (Abeam pic, UK) at 1/200 dilution or with 
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Table 1. Yeast strains and plasmids used in this study. 





Yeast strain or Plasmid 


Genotype 


Source 


BY4743pam 16J::kanMX6/PAM16 


leu2A/leu2A his3A/his3A ura3A/ura3A lys2A/LYS2 metl5A/MET15 
pam 16A::kanMX6/PAMW 


Winzeler, E.A. [15] 


YPH499pam 16J::ADE2-pPAM16 


leu2A ade2A ura3A hls3A trplA lys2A paml6A::ADE2 [pPAM16 URA3] 


Guiard, B. (kind gift) 


BY-pam 75^MAGwt 


leu2A/leu2A his3A/his3A ura3A/ura3A lys2A/LYS2 metl5A/MET15 
pam 16A::kanMX6[pRS41 5-WT-MAGMAS] 


This study 


BY-pam /fi^-IVlAGAs„76Asp 


leu2A/leu2A his3A/his3A ura3A/ura3A lys2A/LYS2 met15A/MET15 
pam 16A::kanMX6 [pRS41 5-Mut-MAGMAS] 


This study 


YPH499-pam /fiz^MAGwr 


leu2A ade2A uraSA his3A trplA \ys2A paml6A::ADE2 [pRS41 5-WT-MAGMAS] 


This study 


YPH499-pam 76z/-MAGas„76asp 


leu2A ade2A ura3A his3A trplA lys2A paml6A::ADE2 [pRS41 5-Mut-MAGMAS] 


This study 


YPH499-pam/&^pPAI\/116-MAGwT 


leu2A ade2A ura3A his3A trplA lys2A paml6A::ADE2 [pPAM16 URA3] [pRS415- 
WT-MAGMASj 


This study 


YPH499-pam J6J-pPAM 1 6-MAGAsn76Asp 


leu2A ade2A ura3A hisSA trplA lys2A paml6A::ADE2 [pPAM16 URA3] [pRS415- 
Mut-MAGMAS] 


This study 


pRS415-WT-MAGMAS 


Amp"CfW LEU2 TEF MAGMAStvr 


D'Silva P. et al. [17], 


pRS4 15-Asn76Asp-MACMAS 


Amp"C£N LEU2 TEF MAGMASas„76A!p 


This study 


pYX232-mtCFP 


Amp"2 micronTPl preSu9 GFP TRPl 


Westermann B. and Neupert 
W. [41] 


pUG34 DsRed.SKL 


Amp''MET25DsRed.SKL HIS3 


Kuravi K.et al. [42] 



doi:l 0.1 371/journal.pgen.l 004311. tool 



Ab specific to Col X (Qiaartett) at dilution 1/20 using DAKO 
En Vision kit. Images were captured with an Olympus PD70- 1X2- 
UCB microscope. All experimental procedures were approved by 
the French Animal Care and Use Committee. 

Plasmids, yeast strains and growth conditions 

Plasmids and yeast strains used in this study are described in 
table 1. The yeast expression vector pRS41 5-WT-MAGMAS, 
encoding full-length human WT (WT) MAGMAS under the TEF 
promoter was provided by D'Silva et al. [17]. The pRS415- 
Asn7 6Asp-MAGMAS plasmid carrying the point mutation p.As- 
n76Asp was generated bom pRS41 .5-WT-MAGMAS thvoTdgh a site 
directed mutagenesis using the QuickChange protocol (Stratagene) 
according to the manufacturer's recommendations. Direct se- 
quencing was performed to confirm the successful introduction of 
the mutation. For construction of yeast strains expressing the 
human WT/mutated MAGMAS, we used Saccharomyces cerevisiae 
Apaml6 strains (provided by Delahodde A). For the yeast 
complementation assay, the diploid strain BY4743 carrying the 
deleted paml6 allele (/Jomiff/liikanMXG) in a heterozygous state 
was first transformed with pRS4 1 5-WT-MAGMAS or pRS415- 
Asn76Asp-MAGMAS plasmids, accordingly and then subjected to 
sporulation and tetrad dissection using standard techniques [39]. 
Spores from a single ascus were first grown at 28°C (301.15°K) in 
YPD media then replica-plated to G4 1 8 plates for the selection of 
spores carrying the deleted paml6 allele and to Leucine drop-out 
plate for the selection of pRS415MAGMAS positive spores. The 
haploids 'SY-paml6A-MAGy^-[ and BY-/)am7ffzl-MAGAs„76Asp 
were thus generated and spotted onto YPD and YPG to analyze 
their phenotype. For the construction of YPH499-/)am76'zl- 

pPAM16-MAGwY and YPH499-/)flm76'zf-/;/yA/76'-MAGA.„76A_sp, 
the haploid YJ'ii499paml 6A ::ADE2 strain containing a plasmid 
encoding the yeast PAM16 gene (kindly provided by Guiard B), 
was transformed either with pRS41 5-WT-MAGMAS or pRS415- 
Asn76Asp-MAGMAS plasmids, respectively then patched onto the 
selection plates (lacking leucine and uracUe). Transformants were 
subsequently rephcated to 5-Fluoroorotic (5-FOA) containing 



plates, recovered and grown at 28°C (301.15°K), thus generating 

YPH499-/;flm76'/l-MAGwT and YPH499-/7flm76'zl-MAGAs„76Asp 
strains that have lost the plasmid containing the yeast PAM 16 
gene and the URA3 marker gene [40] . Ten fold serial dilutions of 
the different generated strains were subsequentiy spotted onto 
YPD plates then incubated at diflFerent temperatures for 3 days. 
Unless otherwise noted, yeast cultures were grown at 28°C 
(301.15°K). YPD (1% bactopeptone, 1% yeast extract, and 2% 
glucose) and YPG (1% bactopeptone, 1% yeast extract, and 2% 
glycerol) were used as rich fermentable and non-fermentable 
growth media respectively. WO (0.67% yeast nitrogen base without 
aa and 2% glucose) supplemented with appropriate nutritional 
requirements according to the strains was used as a minimal 
medium. For solid media, 2% bacto agar (Difco) was added. 

Protein import and western blotting 

For the import activity analysis, Y'P}A4:'d9paml6A-PAM16, YPH- 
paml6A-y^AG^J^l-[ and YPH-/)ami6/f-MAGAsn76Asp cells were first 
grown to ODgoo nm of 1 .0 at 28°C (30 1 . 1 5°K) in rich medium then 
shifted to 37°C (310.15°K), for 8h (28800 s) to induce non- 
permissive conditions. Cell pellets were suspended in ice-cold TCA 
buffer (TrisHCl pH 8 20 mM, Ammonium Acetate 50 mM, 
EDTA 2 mM, protease inhibitor solution, and TCA10%) then 
lysed by adding glass beads and shaking vigorously for 300 s (5 fold 
60 s in ice between). After centrifugation of 600 s (13,000 rpm) at 
4°C (277.15°K), the proteins were solubUized in TCA-LaemmU 
loading buffer (SDS/glycerol stock solution 48%, Tris EDTA 
40%, pi-mercaptoethanol 0.5%, PMSF 0.2%, protease inhibitor 
solution 0.2%) and separated by SDS-PAGE. The proteins were 
subsequently transferred to a nitrocellulose membrane that was 
first blocked, then incubated 1 h at room temperature with a 
blocking solution containing the monoclonal HSP60 Ab (provided 
by Dujardin G.) at 1/15000 dilution. After several washes with 
TBS IX (Tris Buffered Saline IX), the membrane was immersed 
for 1 h at room temperature with HRP - conjugated secondary 
anti-mouse Ab (at 2/15000 dilution in a blocking solution) then 
washed again with TBS IX. Secondary antibodies bound to the 
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membrane were detected by chemiluminescence using the 
Thermo Scientific PierceECL 2 according to the manufacturer's 
instructions. The same protocol was used for the immunodetection 
of MAGMAS protein in YPH499, YFH-pam] 6A-^IAGwr and 
YPH-^amiffzf -MAGAsnyeAspStrains. For this purpose, we incubated 
the membrane with a blocking solution containing the polyclonal 
anti-MAGMAS Ab (Abeam pic, UK) at a 1/1000 fmal 
concentration. Therefore, the secondary Ab used in this case 
was HRP - conjugated secondary anti-rabbit Ab (at 2/15000 
dilution in a blocking solution). Phosphoglycerate kinase (PGK) 
and Porin were detected using anti-PGK (Invitrogen) and anti- 
Porin (Invitrogen), accordingly. 

Fluorescence microscopy 

In order to label mitochondria, YPH-jdamiffzf-MACrvvT and 
YPH-/)am76zl-MAGAsn76Asp strains were transformed with the 
pTX232-mtGFP plasmid that encodes the mitochondrial matrix- 
targeted GFP [41] and grown for 2 days at 28°C (301.15°K). CeUs 
were secondly incubated overnight at the non-permissive temper- 
ature to induce the phenotype. For the visualization of peroxi- 
somes. BY-/)am76zl-MAGv\ I and BY-Zmmiffzl-MAGAsnyBAsp were 
transformed with pUG34DsRi:d.SKL encoding the red fluorescent 
protein DsRed containing the C-terminal peroxisomal targeting 
signal type 1 (PTSl, DsRed-SKL). Transformants were grown for 
2 days at 28°C (301.15°K) in minimal medium. Harvested cells 
were washed in Phosphate-buffered saline IX (PBS IX) then 
transferred to microscopic slides for imaging. For observations at 
high temperature, cells grown 2 days at 28°C (301.1.5°K) in 
minimal medium, were shifted to 37°C (310.15°K), for 6 hours 
then harvested, washed and transferred to slides for visualization. 
The sUdes were examined with a DMIRE2 microscope (Leica, 
Deerfield, IL). For GFP fluorescence, filters used were 450/ 
490 nm for excitation and 500/550 nm for emission whereas for 
Ds.Red fluorescence, filters were 542/582 nm for excitation and 
604/644 nm for emission. Images were captured using a CCD 
camera (Roper Scientific, Tucson, AZ) and Metamorph software 
(Universal Imaging, WestChester, PA) was used to deconvolute 
series and treat the images. 

Web resources 

The URLs for data presented herein are as follows: 
Online MendeHan Inheritance in Man (OMIM), http://www. 
omim.org 

Polyphen, http:/ /genetics.bwh.harvard.edu/pph2/ 
MutationTaster http:/ /www.mutationtaster.org/ 
SIFT http:/ /sift.bii.a-star.edu.sg/ 
PROVEAN http:/ /provean.jcvi.org/index.php/ 
PLINK http://pngu.mgh.harvard.edu/purceU/plink/ 
Primer3 software http://frodo.wi.mit.edu/ 

Supporting Information 

Figure SI Pedigrees and haplotypes of families Fl and F2. 
Markers are reported from telomere (fc/i) to centromere {bottom) on 
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chromosome 16pl3.3. Blackened symbols represent affected individ- 
uals. Disease-bearing chromosome is represented in black. The red 
box represents the minimal ancestral homozygous haplotype 
spanning 1.9 Mb between markers D16S758 and D16S243 and 
shared between F1-IV.3 and F2-IV.3. The normal (A) or mutated (G) 
allele oi MAGMAS \s, indictated within the haplotype. 
(TIF) 

Figure S2 Protein import pathways into mitochondria. (A) Several 
sophisticated transport machin(;ri<;s mi;diate recognition, import, 
sorting and assembly of preproteins into a specific subcompartment 
of the mitochondria. Preproteins targeted to the matrix cross the 
outer mitochondrial membrane through the TOM complex then 
pass through the inner membrane via the TIM23^''^^'^''^. Once 
processed into the matrix, preproteins are folded to their active forms 
after the cleavage of their presequences by the mitochondrial 
processing peptidase (MPP). HSP60 assists in the correct folding and 
assembly of the imported proteins. (B) TIM23'^°™^: the translo- 
case involved in preproteins import into the mitochondrial matrix. It 
is composed of TIM23 translocase associated with the PAM complex 
(In red). The ATP-dependant core of PAM, mtHSP70, drives the 
translocation and the unfolding of preproteins. The regulation of the 
position and activity of mtHSP70 is mediated by five other subunits 
of the PAM complex, consisting of a soluble nucleotide exchange 
factor MGEl and four membrane bound co-chaperones TIM44, 
DNAJC19, PAM16 and PAM17. Via its C-terminal J-like domain, 
PAM 16 (red arrow) interacts with DNAJC19 thus enabling its 
tethering to the translocon. IM, inner membrane; IMS, intermem- 
brane space; OM: outer membrcine. 
(TIF) 

Figure S3 Primers used for expression analysis by quantitative 
Real-Time RT-PCR. F: Forward primer. R: Reverse primer. Fl, 

Rl are specific to the C0R07 transcript; F2, R2 to the MAGMAS 
transcript and F3, R3 to the C0R07-MAGMAS transcript. 
(TIF) 

Table SI The list of common exonic variations shared at a 
homozygous state between both probands F1-IV.3 and F2-IV.3. 
(DOCX) 
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